EVALUATING RADIOMETRIC SENSITIVITY OF LANDSAT 8 OVER 

COASTAL/INLAND WATERS 


Nima Pahlevan 1,2 * , Jianwei Wei 2 , Crystal B. Schaaf, John R. Schott 

'NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
2 Sigma Space Corporation, Lanham, MD 20706, USA 
nima.pahlevan@nasa.gov 

'J 

School for the Environment, University of Massachusetts at Boston, MA 02125, USA 
Jianwei.Wei@umb.edu , crystal.schaaf@umb.edu 

4 Center for Imaging Science, Rochester Institute of Technology, Rochester, NY, 14623, USA 

schott@cis.rit.edu 


ABSTRACT 

The operational Land Imager (OLI) aboard Landsat 8 was 
launched in February 2013 to continue the Landsat’s 
mission of monitoring earth resources at relatively high 
spatial resolution. Compared to Landsat heritage sensors, 
OLI has an additional 443 -nm band (termed coastal/aerosol 
(CA) band), which extends its potential for 
mapping/monitoring water quality in coastal/inland waters. 
In addition, OLI’s pushbroom design allows for longer 
integration time and, as a result, higher signal-to-noise ratio 
(SNR). Using a series of radiative transfer simulations, we 
provide insights into the radiometric sensitivity of OLI when 
studying coastal/inland waters. This will address how the 
changes in water constituents manifest at top-of-atmo sphere 
(TOA) and whether the changes are resolvable at TOA 
(focal plane) relative to OLI’s overall noise 1 . 

Index Terms — OLI, Landsat 8, radiometric 
sensitivity, water constituents, water quality 

1. INTRODUCTION 

The Landsat program has provided the longest record of 
Earth Observation (EO) datasets amongst all of the existing 
EO programs. This valuable asset (which is freely released 
to the public) allows for regular monitoring of changes in 
the earth natural resources over the past four decades. The 
operational Land Imager (OLI) aboard Landsat 8 was 


1 This project was awarded (No. G12PC00072) by the United States 

Geological Survey (USGS) under the contract No. S201 30000023292. 


launched in February 2013 to continue the Landsat’s 
mission for monitoring earth resources at relatively high 
spatial resolution. The enhancements include the addition of 
a 443 -nm channel and a 12-bit radiometric resolution. 

The Landsat datasets, prior to the launch of Landsat 8, have 
been utilized in many studies for mapping water 
constituents. A few studies have demonstrated the potential 
of OLI in mapping optically active components (OACs) of 
water, including the near-surface concentrations of 
chlorophyll-a (CHL) and total suspended solids (TSS), as 
well as the absorption of the colored dissolved organic 
matter (CDOM), i.e., a CD0 M [T 2]. Fig. 1 illustrates the 
scene-derived signal-to-noise-ratio (SNR) curve obtained 
from ~ 30 recorded OLI scenes collected across different 
latitudes (solar zenith angle (0 S ) ranging within 20° — 60°). 
The relatively high SNRs (for the CA and blue bands) 
comparable to those of coastal-ocean missions indicate 
promise for coastal/inland water studies. 

Recognizing the potential of OLI for studying bio -optical 
properties in coastal/inland waters, we examine spectral- 
radiometric (spectroradiometric) sensitivity of OLI for 
studying coastal/inland waters, i.e., how well changes in 
CHL, TSS, and a CD0 M are detectable in the OLI’s spectral 
bands, and how well the changes are resolvable at TOA 
(focal plane) relative to the system noise. 

2. APPROACH 

2.1. METRIC 

Here, we introduce a realistic metric to measure the OLI’s 
spectroradiometric (spectral-radiometric) sensitivity to 



Fig. 1. The average scene-derived SNR (n = 30) calculated over 
uniform bodies of water. The dashed/dashed-dot curves indicate 
maximum and minimum SNRs derived from the scenes. 

changes in water constituents in coastal/inland waters. The 
question to be answered is how well changes in OACs 
(ACHL, ATSS, Aa CD0M ) can be discerned at TO A using the 
OLI instrument. From a radiometric standpoint, it is 
recognized that when the changes in the signal at TO A is 
less than the system noise level, the retrieved near-surface 
properties, e.g., water-leaving radiance; L w , will bear 
significant uncertainties. 

Here, we evaluate the (desired) signal versus the overall 
system noise for each band. The desired signal, is 
represented by the transmitted water-leaving radiance, i.e., 
tL w , where t is the diffuse attenuation for the sea-to-sensor 
path. The desired change in the signal (given the changes in 
CHL, TSS, and 3 CD0 m) manifests in tAL w . From this term, 
the differential signal-to-noise ratio (A SNR) for the OLFs 
effective responses is defined as following [3]: 

ASNR(A 0 ; AOAC) = = 1 w n w 1 (1) 

where N is the effective system noise and represents the 
system overall noise due to photons, striping, quantization, 
electronics, etc. The ASNR metric, in fact, measures the 
signal contrast versus the overall system noise and is 
directly related to the noise equivalent change in radiance 
(NEAL). A 0 denotes OLFs center wavelengths within the 
visible-near-infrared (VNIR) region at 443 nm , 482 nm , 562 
nm , 655 nm , which represent CA, blue, green, and red 
bands, respectively. The differences in the water 
constituents (AOAC) are modeled via Hydrolight [4] 
simulations for two adjacent pixels (p x and p 2 ) in the scene 
located at an arbitrary location. The AOAC parameter 
denotes a change in one of the constituents, i.e., ACHL, 
ATSS, or Aa CDC) M(440). When a desired change in one 
component is desired, the two other components are held 
fixed. For example, a change in CHL (ACHL) is estimated 


when TSS and a CDO M(440) as well as the spectral shapes of 
the inherent optical properties (IOPs) remain constant [4]. 
The band-specific spectroradiometric sensitivity (80 AC) of 
OLI is then estimated at ASNR = 1, i.e., when tAL w = N in 
Eq. 1. In essence, this defines a threshold (detection limit) 
below which detecting changes is limited by the 
instrument’s spectral-radiometric capability. The detection 
limits are referred by6CHL, 6TSS, and Sa CDO M(440), 
respectively. The terms in the numerator (Eq. 1) are 
simulated through radiative transfer modeling while the 
noise term in the denominator is derived from measurements 
of selected OLI scenes. 

2.2. SIMULATIONS 

In this study, water-leaving radiance (L w ) for various water 
types is modeled with the Hydrolight code, whereas the 
effective noise is estimated through scene analysis. 

The Hydrolight simulations were conducted for various 
CHL, TSS, and a CDO M(440) levels within the 0.2 to 10.0 
mg/m 3 , 0.0 to 12.0 g/m 3 , and 0.05 to 2.0 m -1 ranges, 
respectively. The input concentrations span a reasonable 
range representing moderately productive, mesotrophic 
waters to highly turbid/trophic waters. The VNIR region 
(400 to 900 nm) in-water simulations at 5 nm resolution 
were conducted for three solar geometries, i.e., 0 S = 
25°, 45°, 60°, and the Julian day of 270 (representing an 
average Earth-Sun distance). The Fournier-Forand phase 
function with 2% backscattering fraction (b b /b = 0.02) 
was chosen for all the simulations . It should be noted that a 
moderate wind speed, i.e., 5 m/s , and aerosol optical 
thickness of 0.1 at 550 nm were used for the simulations. 
While it is desired to perform simulations with very fine 
OAC increments, a compromise between the computation 
time and the simulation step sizes was made: CHL was 
chosen to have increments of 0.2 units of concentrations, 
TSS was incremented by 0.1 units, and the reference CDOM 
absorption, i.e., a CD0M (440) = 0.18 m _1 , was scaled by 
multiplicative factors ranging from 0.3 to 11 with 
increments of 0.06. This sampling scheme provides 
sufficiently fine-scale samples for properly analyzing OLFs 
spectroradiometric sensitivity. The Hydrolight-modeled 
water-leaving radiance (L™) were then propagated to the 
TO A using MODTRAN [5] to obtain simulated TO A 
radiance (L™) and tL™. The TOA quantities were then 
passed through the OLFs RSR functions. Note that the 
relative spatial response (RSR) functions represent the 
cascade of the performance of all system components (e.g. 
quantization, optics, and electronics). 

2.3. SCENE-DEPENDENT NOISE 

To estimate the effective scene noise (N) on the focal plane, 
the scene-derived standard deviations and mean TOA (L t ) 
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Fig. 2. The scene-derived noise plotted against TO A radiance for 
the CA band. The linear relationship provides a realistic estimate of 
system noise for a given TOA radiance (derived from OLI scenes). 

radiances were calculated over spatially uniform bodies of 
waters, i.e., 10 X 10 pixels. The relationships constructed 
between N and L t for OLI scenes (n = 30) are used to 
deduce the noise level (N) for the model-derived TOA 
radiances (L™). Fig. 2 illustrates the observed signal-noise 
relationship for the CA bands. The diversity of the selected 
scenes for noise and signal estimation provides an extended 
dynamic range of observed total sensor-reaching radiances. 
This allows for evaluating the effective noise for various 
modeled TOA radiances. Linear regression analysis is 
applied to the modeled TOA radiances to estimate 
equivalent effective noise, and, as a result, ASNR (Eq. 1). 

3. RESULTS 

Many combinations of OACs were examined to map out 
OLI’s spectroradiometric sensitivity to changes in OACs for 
a nominal solar zenith angle, i.e., 6 S = 40°. Later, the 
impact of changes in 6 S will be discussed. Fig. 3 illustrates 
the band-specific CHL detection limits (SCHL) for three 
different ranges of CHL, i.e., 0.1 < CHL < 1, 1 < CHL < 
5, and 5 < CHL < 12, referred to as case I, II, and III, 
respectively (not to confuse with widely used case I and 
case II water types). Realistic TSS and a CDO ]vi(440) values 
in accordance with each CHL range were incorporated. It is 
revealed that for all cases, the blue band (482 nm) shows the 
highest sensitivity to changes in CHL, i.e., on average, 
SCHL(482) > 0.5 mg/m 3 . This implies difficulties in 
detecting changes smaller than 0.5 units CHL on the focal 
plane using this single band. The spectroradiometric 
sensitivity of the blue band degrades for the most trophic 
waters studied here (5 < CHL < 12). This applies to all 
OLI’s bands. While the CA and the green bands, on 
average, exhibit similar sensitivity to the changes in CHL, 
the CA band is slightly better for case I (least trophic here) 
and case II (moderately trophic) waters. Nevertheless, in 
highly trophic waters, the green band shows lower detection 



Fig. 3. The OLI’s band-specific detection limits for CHL. The CA 
band (443 nm) is found as sensitive as the green band for low CHL 
conditions. Three CHL level are as following 0.1 < CHL < 1, 

1 < CHL < 5, and 5 < CHL < 12. Note that the lower 6CHL the 
more sensitive is OLI (relative changes in CHL). 

limits. In relatively low-productivity waters, the red band is 
nearly as sensitive as the blue band (8CHL( 655)~ 0.5), 
whereas it represents similar detection limits as in the green 
band for more turbid waters (case II and III). 

The spectroradiometric sensitivity of the OLI bands with 
respect to the changes in TSS is shown in Fig. 4. The water 
types are classified into three broad regions, 0.1 < TSS < 1, 
1 < TSS < 5, and 5 < TSS < 12. Compared to CHL, OLFs 
sensitivity relative to ATSS is more than 10 x better than 
that for ACHL. As expected, the sensitivity diminishes, 
particularly in the CA band, for highly turbid waters (with 
high CDOM). The green band, on average, shows most 
sensitivity to changes in TSS. This is attributed to the 
relatively small total absorption within this region of the 
spectrum. Overall, the OLI sensor is expected to be capable 
of detecting differences in TSS larger than, on average, 
0.1 g/m 3 , i.e., STSS « 0.1, in different water types. The 
water types are classified into three broad regions: 0.1 < 
TSS < 1, 1 < TSS < 5, and 5 < TSS < 12. Compared to 
the sensitivity to CHL, the sensitivity with respect to TSS is 
found to be 10 x better. As expected, the sensitivity 
diminishes, particularly in the CA band, for highly turbid 
waters (with high CDOM). The green band, on average, 
shows most sensitivity to changes in TSS. This is attributed 
to the relatively small total absorption within this region of 
the spectrum. 

The detection limits associated with a CDO M(440) for the 
OLFs CA, blue, and green bands are illustrated in Fig. 5. 
The red band was found insensitive to the changes in 
a CDOM(440). In a similar fashion as for CHL and TSS, the 
water types were classified into three broad regions, 
0.01 < acDOM(440) < 0.1, 0.1 < a CDO M(440) < 0.5, and 
0.5 < a CDC )M(440) <2. In relatively low CDOM 
concentrations, i.e., a CDO ]vi(440) < 0.5 m -1 , with CHL and 
TSS remaining below 5 units, the blue and the green bands 
exhibit the highest sensitivity (6a CDC >M(440) < 0.1), 



Fig. 4. The OLI’s band-specific detection limits for TSS. The 
sensor is, on average, sensitive to ATSS > 0.1 g/m 3 in nearly all 
bands and for all water types. The CA band is found as sensitive 
as the other bands in slightly turbid waters.. Three TSS level are as 
following 0.1 < TSS < 1, 1 < TSS < 5, and 5 < TSS <12. 


Fig. 5. The OLI’s band-specific detection limits for a CD0M ( 440). 
The sensor is, on average, sensitive to a CD0M ( 440) >0.1 m -1 
for water types with a CD0M ( 440) < 0.5 m -1 . The red band is 
excluded from this chart due to its insensitivity to changes in 
CDOM. Three CDOM levels are as 0.01 < a CD0M ( 440) < 0.1, 
0.1 < &cdom( 440) < 0.5, and 0.5 < clcdom^ 40) < 2. 


whereas the sensitivity of the CA band is just above 
0.15 m -1 . The detection limit for CDOM-rich waters, i.e., 
0.5 < a CDC) M(440) < 2, deteriorates for the OLFs CA and 
blue bands (due to significant CDOM-induced absorption in 
the ultraviolet-blue portion of the spectrum). In such waters, 
the green band provides the best sensitivity to changes in 
a CDOM(440), i.e., 6a CDO M(440)~0.1. One must bear in 

mind that the green band is also relatively highly sensitive to 
TSS and CHL and care must be taken when evaluating the 
responses in this band. 

The detection limits can be improved n x by averaging the 
signal over n x n areas. The average detection limits of 0.5 
units for CHL, i.e., SCHL = 0.5 mg/m 3 , for instance, will 
be enhanced to SCHL = 0.167 mg/m 3 , if a 3x3 
smoothing kernel is applied on the recorded OLI datasets. 
One must balance the loss of spatial details against the gain 
in the retrieval sensitivity in areas where spatial variability 
in OACs is deemed minimal, i.e., smoothing will improve 
retrieval performance. For high-latitude OLI observations 
(not presented here), where 6 S « 60°, our analysis indicates 
that there is, on average, a 30% reduction in the detection 
limits; on the other hand, when 9 S = 25°, the radiometric 
sensitivity is, on average, only improved ~6%, when 
compared to the mid-latitude detection limits (Figs. 3,4, and 
5). Note that the A SNR metric is intended as an engineering 
tool to assess the expected performance limits. Actual 
retrievals will be improved by the use of multiple bands and 
may be degraded by calibration errors and less ideal 
atmospheric compensation technique. 

4. CONCLUSION 

The analysis of the spectroradiometric sensitivity of OLI 
using A SNR as metric gives clues on the errors anticipated 
for the retrieval of the biogeochemical properties (CHL, 


TSS, a CD0 ]yi)- F or an average solar zenith angle (0 S = 45°), 
it was found that OLI is, on average, sensitive to changes in 
CHL, TSS, and a CDO M(440) larger than 0.5 mg/m 3 and 
0.1 g/m 3 , and 0.1 m -1 , respectively. For high-latitude 
areas with 0 S ~ 60°, the sensitivity degrades up to 30%. To 
improve the retrieval of OACs over large bodies of waters, 
the OLI scenes may be smoothed to lower the noise 
contribution. This study demonstrated OLFs band-specific 
sensitivity to changes in OACs, despite its relatively broad 
spectral bands (when compared to the corresponding ocean 
channels of MODIS and those of VIIRS). In addition to 
limitations in the OLFs radiometric sensitivity, further 
complications are introduced when correcting for the 
atmospheric effects. 
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